Minimum uncertainty measurements of angle and angular momentum 
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The uncertainty relations for angle and angular momentum are revisited. We use the exponential 
of the angle instead of the angle itself and adopt dispersion as a natural measure of resolution. We 
find states that minimize the uncertainty product under the constraint of a given uncertainty in 
angle or in angular momentum. These states are described in terms of Mathieu wave functions and 
may be approximated by a von Mises distribution, which is the closest analogous of the Gaussian 
on the unit circle. We report experimental results using beam optics that confirm our predictions. 
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Light carries and transfers energy as well as linear and 
angular momentum. The angular momentum contains a 
spin contribution, associated with polarization, and an 
orbital component, linked with the spatial profile of the 
light intensity and phase 0]. 

The seminal paper of Allen et al 0] firmly establishes 
that the Laguerre-Gauss modes, typical of cylindrical 
symmetry, carry a well-defined angular momentum per 
photon. In the paraxial limit, this orbital component 
is polarization independent and arises solely from the az- 
imuthal phase dependence e 4m< ^, which gives rise to spiral 
wave fronts. The index m takes only integer values and 
can be seen as the eigenvalue of the orbital angular mo- 
mentum operator. In consequence, the Laguerre-Gauss 
modes constitute a complete set and can be used to rep- 
resent quantum photon states 0, 0, • 

The possibility of exploiting these light fields for driv- 
ing micromachines, and their applications as optical 
tweezers and traps, have attracted a good deal of at- 
tention |E S ■ Moreover, entangled photons prepared 
in a superposition of states bearing a well-defined or- 
bital angular momentum provide access to multidimen- 
sional entanglement. This is of considerable importance 
in quantum information and cryptography because with 
these states more information can be stored and there is 
less sensitivity to decoherence 0, 0, 0, 0, . 

Recently, by precise measurements on a light beam, an 
experimental test of the uncertainty principle for angle 
and angular momentum has been demonstrated 0, ^3 • 
The idea is to pass the beam through an angular aper- 
ture and measure the resulting distribution of angular- 
momentum states Moreover, one can even identify 
the form of the aperture that corresponds to the mini- 
mum uncertainty states for these variables. 

In the following, we deal with cylindrical symmetry; 
we are concerned with the planar rotations by an angle 
(f> in the plane x-y, generated by the angular momentum 
along the z axis, which for simplicity will be denoted 
henceforth as L. In this respect, we recall that the proper 
definition of angular variables in quantum mechanics is 
beset by well-known difficulties [171 llSj . For the case of 



a harmonic oscillator, the problems essentially arise from 
two basic sources: the periodicity and the semibound- 
edness of the energy spectrum. The first prevents the 
existence of a phase operator in the infinite-dimensional 
Hilbert space, but not of its exponential. The second 
entails that this exponential is not unitary. 

Although we have here the same kind of problems 
linked with the periodicity, the angular momentum has a 
spectrum that includes both positive and negative inte- 
gers, which allows us to introduce a well behaved expo- 
nential of the angle operator, denoted by E 0. Since 
the angle is canonically conjugate to L, we start from the 
commutation relation [isl | 



[E,L]=E. 



(1) 



The goal of this Letter is precisely to develop a com- 
prehensive approach to the minimum uncertainty states 
associated with the relation (JJJ . Our results will corrob- 
orate that the use of E provides a good description of 
the angular behavior and the associated minimum un- 
certainty states turn out to be Mathieu beams in wave 
optics. This will establish a proper basis for information 
processing with these conjugate variables. 

We also stress that, since angle is 2n periodic, the cor- 
responding quantum statistical description should also 
preserve this periodicity. Provided that a non-periodic 
measure of the angular spread is used, as for example 
the variance, such a resolution depends on the 2ir win- 
dow chosen. To prevent this, we recall that another ap- 
propriate and meaningful measure of angular spread is 
the dispersion [2fj| 
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and P(4>) is the angle distribution. As expected, it pos- 
sesses all the good properties: it is periodic, the shifted 
distributions P(<f> + </)') are characterized by the same 
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resolution, and for sharp angle distributions it coincides 
with the standard variance since |(e^)| 2 ~ 1 + (<f> 2 )- In 
consequence, the statistics of the exponential of the angle 
provides a sensible measure of the angle resolution. 

The action of the unitary operator E in the angular 
momentum basis is 



Elm) = \m-l) 



(4) 



we can restrict ourselves to states with (L) = 0, since 
we readily obtain solutions with mean angular momen- 
tum m by multiplying the wave function by exp(ifh(j)) . 
Alternatively, we observe that the change of variables 
exp(i/i0)4'((/)) eliminates the linear term from @. In ad- 
dition, we can take q to be a real number, since this only 
introduces an unessential global phase shift. We there- 
fore look at solutions of 



where the integer m runs from — oo to +oo. Therefore, 
E possesses a simple optical implementation by means of 
phase mask removing a charge + 1 from a vortex beam. 
The normalized eigenvectors of E are 
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and, in the representation they generate, we can write 
d 



E = 



(6) 



which formally verify the fundamental relation QJ . 

Let us turn to the corresponding uncertainty relation. 
When the standard form (Ai) 2 (AB) 2 > |([i,B])| 2 /4 is 
applied to Eq. Q and the previous notion of dispersion 
is used, we get 



D 2 (AL) 2 > ^(1 -£> 2 ) 



(7) 



This can be recast in terms of the cosine and sine oper- 
ators, C = (E + £t)/2 and S = (E - E^)/2i, yielding 



(AC) 2 (AL) 2 >i|(C)| 2 , 



(A5) 2 (AL) 2 >i|<5)| 2 . 



(8) 

States satisfying the equality in an uncertainty relation 
are sometimes referred to as intelligent states 21]. How- 
ever, in the case of Eq. (JJJ, the inequality cannot be sat- 
urated (except for some trivial cases), since this would 
imply to saturate both relations in JSJ simultaneously. 
In other words, the formulation J7J) is true but too weak. 

To get a saturable lower bound we look instead at 
normalized states that minimize the uncertainty prod- 
uct D 2 (AL) 2 either for a given D 2 or for a given 
(AL) 2 . These have been called constrained minimum 
uncertainty-product states [l4| . We approach this prob- 
lem by the method of undetermined multipliers. The 
linear combination of variations [whether we minimize 
D 2 (AL) 2 for a fixed D 2 or for a fixed (AL) 2 } lead to the 
basic equation 



[L 2 + [lL + (q*E + qE^)/2\YS>) — a\^f), 



(9) 



where /i, q, and a are Lagrange multipliers. We shall 
solve this eigenvalue equation in the angle representation 
fy(tf>) = (4>\*&)- Note first that, without loss of generality, 



d 2 *(</») 



(a - qcos(f>) = 0. 



(10) 



To properly interpret this eigenvalue problem we intro- 
duce the rescaled angular variable 77 = <f>/2, so that 



d 2 ^(i]) 
drj 2 



2qcos(2rj)\ #(77) = 0, 



(11) 



which is the standard form of the Mathieu equation [22( • 
The variable ij has a domain < 77 < 2ir and plays the 
role of polar angle in elliptic coordinates. In our case, 
the required periodicity imposes that the only accept- 
able Mathieu functions are those periodic with period of 
7r or 2ir. The values of a in Eq. (|llfl that satisfy this 
condition are the eigenvalues of this equation. We have 
then two families of independent solutions, namely the 
even and the odd angular Mathieu functions: ce n (r),q) 
and se„(77,(7) with n = 0,1,2,...., which are usually 
known as the elliptic cosine and sine, respectively. For 
ce n (r),q) the eigenvalues are denoted as a n (q), whereas 
for se n (rj,q) they are represented as b n (q): they form an 
infinite set of countable real values that have the property 
ciq < 61 < ai < 62 < • ■ ■■ The parity and periodicity of 
these functions are exactly the same as their trigonomet- 
ric counterparts; that is, ce n (rj,q) is even and se n (r],q) 
is odd in 77, and they have period tt when n is even or 
period 2tt when n is odd. 

Since the 2n periodicity in <f> requires 7r periodicity in 
77, the acceptable solutions for our eigenvalue problem are 
the independent Mathieu functions of even order 



1 2 I ce 2n(v,Q), 

fcanfa.ff) = \j- 

1 se 2n (r],q), 



n = 0,l,. 



(12) 



where the numerical factor ensures a proper normal- 
ization, according to the properties of these functions. 
In what follows we shall consider only even solutions 
ce2n(?7, ?)? although a parallel treatment can be done for 
the odd ones. After some calculations, we get 



(AL) 2 2n = i£dr, 



drj 



ce 2n (?7,a) 



1 



[A 2n (q)-2qe 2n (q)], 



(13) 
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dv ce 2 2n (r],q) cos(2?7) 
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where we have expanded the periodic functions ce2n(?y, q) 
in Fourier series 

oo 

ce 2 „(r?, q)=J2 (?) cos(2fc77), (14) 

fc=0 

and we have integrated term by term, in such a way that 



&2n(q) = A 



(2n) 



(<z)4 2Tl) (?)- 



The coefficients A^^ determine the Fourier spectrum 
and satisfy recurrence relations that are easily obtained 
by substituting l|14(l in the Mathieu equation and can be 
efficiently computed by a variety of methods 23] . 

If we expand ce2n(f?, q) in powers of q and retain only 
linear terms |22j . we have 



■J2Ag n \ q )Ag%(q). (15) 

fc=0 



,a trt ( 2 ") 2 4n 4 -3n 2 + l 2 
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(16) 



4(4n 2 - l) 2 



which shows a quadratic increasing with q of the angular- 
momentum variance and a decreasing of the angle disper- 
sion. The uncertainty product, up to terms q 2 , reads as 



D\ n (AL)l n = n 2 + i [(4n 4 - 5n 2 + 1)(1 - DfJ]. (17) 

It is clear that this product attains its minimum value for 
the fundamental mode n — 0, which saturates the bound 
in Eq. Q for this range of values of q. 

Note that for large dispersions (q — > 0) the fundamen- 
tal wave function may be approximated by 



P {(j)) tx |ce (?7,g)| 2 ~ exp(-qcos< 



(9-0), (18) 



which is the von Mises distribution, also known as the 
normal distribution on the unit circle |2^j . This remark- 
able result shows that optimal states are very close to 
Gaussians on the unit circle j2f|. Curiously enough, it 
has been recently found that the von Mises distribution 
maximizes the entropy for a fixed value of the disper- 
sion HI]. 

In the opposite limit of small dispersions (q — > oo), one 
can also check that 



P {4>) oc \ce (r),q)\ 2 ~exp(- v / gcos< 



(<? 



(19) 

Therefore, von Mises wave functions constitute an excel- 
lent approximation to the minimum-uncertainty Mathieu 
wave functions, except perhaps for intermediate values of 
the dispersions, where a deviation may occur. In Fig. ^ 
we have plotted D (AL) in terms of D: the solid line 
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FIG. 1: Theoretical and experimentally observed uncertainty 
products as a function of the dispersion. Solid line represents 
the fundamental Mathieu beam while the broken line repre- 
sents the von Mises approximation. The difference between 
these two wave functions appears plotted in the inset. We 
have included also the ideal bound given by Eq. 0. 



represents the fundamental Mathieu beam, which pro- 
vides the optimal angular resolution, while the broken 
line represents the von Mises approximation. The very 
small difference between these two curves is magnified 
in the inset. For the purposes of comparison, the ideal 
bound coming from Eq. (JJJ is plotted as a dotted line. 
We stres that the minimum uncertainty states with large 
dispersions present wide angular distributions and vice 
versa. 

This theoretical approach can be experimentally re- 
alized, although our capabilities to prepare states and 
perform measurements on demand are limited by the 
present technology. Figure [21 shows our experimental 
setup. Two spatial light modulators (SLM) were used: 
the amplitude SLM (CRL Opto, 1024 x 768 pixels) gener- 
ates the angular-restricted light beams, while the phase 
SLM (Boulder, 512 x 512 pixels) works as an analyz- 
ing hologram. The beam generated by an Ar laser (514 
nm, 200 mW) is spatially filtered, expanded and colli- 
mated by the lens Li and impinges on the hologram 
generated by the amplitude SLM. The bitmap of the 
hologram is computed as an interference pattern of the 
signal beam U s and an inclined reference plane wave 
Up = Up exp[— ik(x sin 7 + ZCOS7)], where k denotes the 
wave number and 7 is the angle of the transversal compo- 
nent of k with respect to the z axis. After illuminating 
the hologram with the collimated beam, which can be 
approximated by the plane wave Ui — Ui exp(—ikz), the 
field behind the SLM can be written as Z7 t = Ui | U s + U p | 2 . 
The Fourier spectrum of this transmitted beam is local- 
ized at the back focal plane of the first Fourier lens FLi 
and consists of three diffraction orders (—1,0,-1-1). The 
undesired and — 1 orders are removed by the spatial fil- 
ter. After inverse Fourier transformation, performed by 
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FIG. 2: Experimental setup for the generation of beams with 
a von Mises distribution and subsequent detection of the as- 
sociated angular-momentum components. 

the second Fourier lens FL2, a collimated beam with the 
required complex amplitude profile U s is obtained. This 
field impinges on the reflecting phase SLM. The hologram 
on the SLM is of the form Uh = exp[im(f) + iS(x)], where 
4> is the azimuthal angle and S(x) = mod (x,A x ) is the 
sawtooth function of period A x , which ensures deflection 
of the undesired orders coming from the pixel structure 
of the SLM. 

When the field impinging on the phase SLM has a 
structure X)^L-oc a « exp(m</>), the hologram modifies its 
hclicity so that the reflected beam can be written as 
Yl^—-^ a n exp\i(n + m)<p]. The components for which 
n + m = present no helicity and their intensity pro- 
file consists of a light spot localized at the center of the 
beam. The other components have a nonzero helicity and 
the center of the intensity pattern remains dark. From 
this fact, the weight coefficients of the superposition can 
be determined by selective intensity measurements per- 
formed by a pinhole and a power meter. 

A Laguerre-Gauss beam was used to align the setup 
and subsequently a beam with a von Mises distribution 
(|18|l of the transversal amplitude was generated. Each 
angular amplitude width was scanned for values of the 
helicities from n = — 20 to n = 20. Experimentally mea- 
sured uncertainty products are depicted in Fig.^by solid 
circles. 

Given the accuracy of the measurements (indicated by 
error bars in Fig. [I}, they fit quite well the theoretical 
predictions. Our present experiment distinguishes be- 
tween the uncertainty product of optimal states and the 
ideal limit, ft is, however, not possible to dicriminatc 
between the Mathicu and von Mises beams. Keeping in 
mind that von Mises and Mathieu states play the same 
role for the spatial degrees of freedom as Gaussian states 
for quadratures, the observation of the nonclassical be- 
havior of angle and angular momentum is a challenging 
problem left for future studies. 

In conclusion, we have formulated the uncertainty rela- 
tions for angle and angular momentum based on disper- 
sion as a correct statistical measure of error. The optimal 
states were derived and identified with the Mathieu wave 
functions. An optical test of the derived uncertainty re- 
lations was proposed and performed experimentally 
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